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Abstract. We propose a trapped ion scheme en route to realize spin Hamiltonians 
on a Kagome lattice which, at low energies, are described by emergent Z 2 gauge fields, 
and support a topological quantum spin liquid ground state. The enabling element in 
our scheme is the hexagonal plaquette spin-spin interactions in a 2D ion crystal. For 
this, the phonon-mode spectrum of the crystal is engineered by standing-wave optical 
potentials or by using Rydberg excited ions, thus generating localized phonon-modes 
around a hexagon of ions selected out of the entire two-dimensional crystal. These 
tailored modes can mediate spin-spin interactions between ion-qubits on a hexagonal 
plaquette when subject to state-dependent optical dipole forces. We discuss how 
these interactions can be employed to emulate a generalized Balents-Fisher-Girvin 
model in minimal instances of one and two plaquettes. This model is an archetypical 
Hamiltonian in which gauge fields are the emergent degrees of freedom on top of the 
classical ground state manifold. Under realistic situations, we show the emergence of 
a discrete Gauss’s law as well as the dynamics of a deconfined charge excitation on a 
gauge-invariant background using the two-plaquettes trapped ions spin-system. The 
proposed scheme in principle allows further scaling in a future trapped ion quantum 
simulator, and we conclude that our work will pave the way towards the simulation of 
emergent gauge theories and quantum spin liquids in trapped ion systems. 
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1. Introduction 

Topological quantum spin liquids are fascinating states of matter supporting topological 
order and exotic excitations with fractional statistics [1]. Recently, there has been 
intense theoretical activity aimed at underpinning their properties, and identifying 
microscopic Hamiltonians which could support such states. This activity is mostly 
focused on the weak and strong insulator scenarios [1]: in the hrst case, the focus 
is on SU(2) invariant spin Hamiltonian on frustrated geometries [2], which could 
potentially be relevant for a series of solid state materials. In the strong insulator 
scenario, one is instead interested in Hamiltonians displaying exotic constrained 
dynamics [1, 3, 4]. This second route has been particularly fruitful at the theory 
level, as the corresponding model Hamiltonians are amenable to powerful analytical 
techniques and, most importantly, provide a very clean interpretation of quantum spin 
liquids in a gauge theory framework [1, 5]. However, due to the exotic form of the 
basic constraints necessary to realize a strong insulator dynamics, many paradigmatic 
Hamiltonians in this class are still thought of as idealized dynamics without any direct 
physical counterpart. 

Here, we show how a minimal instance of an archetypical example of spin 
Hamiltonian supporting a topological quantum spin liquid ground state, the Balents- 
Fisher-Girvin (BFG) model [6], can be realized in 2D laboratory trapped ion systems. 
The key element is to engineer constrained dynamics in a controlled fashion, and in 
particular, to realize exotic ’plaquette’ constraints typical of spin Hamiltonians in the 
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strong insulator scenario. The BFG model is defined on a Kagome lattice, where each 
site hosts a spin-1/2 degree of freedom, and is 

// = jJ5^sA + T ^ + ftx.). (1) 

VieO / {ij) 

Here, S'/ (a = -|-, —, z) are spin operators acting on the site j, the first term with J^(> 0) 
represents Ising interactions around each hexagonal-plaquette, and the second term 
(with J±) describes the nearest-neighbour spin exchange interaction (we consider here 
the variant of the BFG discussed in Ref. [7, 8], where only nearest-neighbor exchange 
is included). The underlying source of frustration in this model is the emergence of 
a macroscopic classical ground-state degeneracy due to the local hexagonal plaquette 
constraints imposed by the first term in the Hamiltonian, i.e., S'* = 0. Hence, 

the ground state configurations are such that each hexagon in a Kagome lattice consists 
of three up spins and three down spins. This is reminiscent of the ice rule in spin- 
ice models [3, 4] in which, close to each vertex of a square lattice, two spins point 
inwards, and two spins outward. Once quantum fluctuations are turned on, J± < Jz, 
the system is effectively described by a quantum dimer model on the dual (triangular) 
lattice [5, 9]. Thus, in contrast to the 2D version of spin Ice, whose dynamics is described 
by a compact U(l) gauge theory, the low-energy physics of the BFG model is different, 
as the underlying lattice on top of which the quantum dimer model is defined is not 
bipartite. In those cases, the emergent degrees of freedom are Z 2 dynamical gauge 
fields [1], which undergo deconfinement and stabilize a topological quantum spin liquid 
state [6, 8, 10]. 

Our proposal to implement the BFG model is based on the remarkable achievements 
in trapped ion technology in the last decades [11, 12]. These highly controllable quantum 
systems have led to a variety of proposals and realizations of quantum simulators, 
such as the topological hexagonal Kitaev model [13], fermionic lattices [14], SU(2) 
Ising models [15, 16, 17, 18, 19] -including frustrated magnetism [20, 21, 22] and, very 
recently, to the observation of entanglement dynamics in spin chains [23, 24]. Besides 
these studies, which are implemented in linear-(ID) ion crystals, quantum simulation of 
various spin models has also been proposed in 2D ion-crystals [25, 26, 27] and various 
experimental approaches are studied, such as 2D Paul traps [28, 29, 30, 31], Penning 
traps [32, 33], micro- [34] and multi zone- [35, 36] trap arrays. Recent work in a Penning 
trap has demonstrated controllable spin-spin interactions between a few hundred ionic 
spins [19], whereas studies in Paul traps show excellent prospects for implementing such 
interactions as well [30]. 

Here, we show how to generate hexagonal-plaquette spin-spin interactions by (anti) 
pinning appropriate ions in a 2D ion crystal by external potentials which we describe 
below. The spin-spin interactions are mediated by a spin dependent optical dipole force 
interacting with a localised phonon mode that appears as a result of this pinning. This 
localised phonon mode involves all ions in the plaquette but rapidly decays outside it. 
The resulting spin-spin interactions on the hexagonal-plaquette represent a key building 
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block for quantum simulators aimed at studying the fundamental nature of frustrated 
quantum magnetism with emergent gauge helds [2, 4, 37] and may open up new 
directions in ion-based quantum simulators. By pinning multiple ions, more individual 
plaquette interactions can be engineered thus allowing for scaling up the quantum 
simulator. In the present paper, we discuss the minimal instance of a two-plaquette 
implementation, which should lie within current experimental reach. Note that our 
current studies also £t well into the recent trends of cold atom-investigations on synthetic 
gauge helds [38], in particular on dynamical gauge helds [39, 40, 41, 42, 43, 44, 45]. 

The ion pinning can be accomplished by a standing wave light held focused onto the 
ion to be pinned. Such optical potentials can reach curvatures corresponding to trap- 
frequencies in the MHz range. The proposed setup resembles those encountered in recent 
studies involving ions interacting with strong light-helds, such as in anomalous dihusion 
of an ion [46], preparing nonclassical motional states [47], and, in particular, modifying 
the trapping potential locally in a rf-trapped ion crystal for quantum simulations [48, 49]. 
Additionally, optical trapping of an ion has been shown experimentally using either a 
single beam dipole trap [50, 51, 52] or an optical lattice [53, 54, 55, 56]. An interesting 
alternative option would be to dress the ion of interest with a Rydberg state, such that 
the dipole moment induced by the ionic trapping held causes a change in local trapping 
frequency for the ion [57, 58]. Here, no spatial variation of the pinning laser would be 
required and additional microwave helds could be employed for hne-tuning [59, 60]. 

The paper is organized as follows: in section 2 we discuss the physical setup of 
the ion-laser system and the governing atom-light interaction Hamiltonian. The spin 
models are derived from the interaction Hamiltonians, we show how to generate spin- 
spin hexagonal interactions in a single plaquette using N = 7 and two plaquettes 
using = 19 ion crystals. In section 3 we discuss the ground state properties and 
magnetization dynamics for the single and double plaquettes spin-systems, and the 
results are compared with that of the original BFG model. Realistic numbers for the 
laser parameters for a calcium ion-setup are given in section 4. Finally, we summarize 
the paper with an outlook in section 5. 

2. The Balents-Fisher-Girvin model in a trapped ion quantum simulator: 
building blocks 

2.1. Spin-spin interactions in a two-dimensional ion crystal 

A sketch of our proposed experiment is shown in hgure 1. The system consists of a self- 
assembled 2D ion crystal confined in the xy plane by what we, for simplicity, assume 
to be a radially symmetric trap with trap-frequencies: ^ ^z- The internal level 

structure of the ions consists of three long-lived low-lying atomic states, these can be 
either three Zeeman states in the ground-state manifold or a combination of metastable 
and ground states. The states |2) =|J,) and |3) =|t) encode the spin-half states of the 
BFG model. The ion crystal supports quantized phonon modes that are used to transmit 
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Figure 1. (a) The schematic picture of the ion-laser system for generating 
hexagonal-plaquette interactions. A pair of counter-propagating Raman fields 
with wave vectors (frequencies) kiz (wi) and k 2 Z (a; 2 ) are used for inducing 
spin-spin couplings. OLl and 0L2 (yellow in colour) are two optical lattices 
focused on ions 1 and 2 to modify their transversal (along the 2 ; axis) trapping 
frequencies locally, (b) The level scheme: the internal level structure consists of 
three long-lived low-lying states and a manifold of excited states. The Raman 
lasers (blue and red arrows), detuned by 6m with respect to the frequency of 
a motional mode, generate state dependent optical shifts in the spin states 
|2) = I I) and |3) = | t) with an energy separation of which simulate the 
spin-spin interactions of the type 5* (8> Si between different ions i and j. The 
ions 1 and 2 prepared in state |1) experience an additional optical lattice that 
is created by very far detuned laser fields (yellow arrow) from the excited state 
|e) with a detuning Aql- The lattices modify their local trapping frequencies 
along the 2 : axis. Additional laser fields, used for generating quantum dynamics 
in the spin system, are not shown in the figure for the sake of clarity. 


spin-spin interactions between the ions via a Raman pair of counter-propagating laser 
fields as shown in hgure lb. Ions sitting in the middle of an hexagonal sublattice (see, 
e.g., ions 1 and 2 in Fig. 1) are prepared in the state |1) and strongly pinned by a 
standing wave laser held that is far detuned from any excited states {|e)} to locally 
modify the phonon spectrum. 

The details of laser mediated spin-spin interactions in trapped ion systems have 
been extensively described elsewhere [11] and are briehy discussed in Appendix A. 
Two counter propagating laser helds are far detuned from the transitions |J,) ^ \e) 
and It) |e). When the Rabi frequencies of these lasers are much smaller than the 
detunings, the excited state |e) can be adiabatically eliminated. To engineer spin- 
spin interactions of the form ® Sz, the frequency difference of the laser beams uj 
is tuned close to the phonon frequencies Um, where m denotes the particular phonon 
mode and 5m = — ^m- Throughout the paper we assume to be in the Lamb-Dicke 

regime, + 1) "C 1, where = Qmobmki are the Lamb-Dicke parameters, Um 
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the number of phonons in mode m, qmo = \/h/2MuJm with M the mass of an ion 
and hm = {&m} f^e (normalised) phonon eigenvectors. As long as rimO,j{n) 6^., 
with the two photon Rabi frequency Qj{n) = Qi{n)Q 2 {'n)/A for each laser and state 
n =t, i and Ai A 2 S> f2i,2(n) the detuning for each laser, the Raman laser only 
transiently excites phonons and we can integrate out the phonon degrees of freedom 
altogether to obtain effective spin-spin interactions between ion i and j of the form 
Hzz = Yhij ^ coupling matrix 

jij ^ ^ ^2) 

m=l 

with fl] = The strength and range of Ej are determined by the Rabi 

frequencies Q], the detuning 6m and the amplitude of oscillation of each ion in the 
mode through the Lamb-Dicke parameters t]1^. Note that we only assumed to be in the 
Lamb-Dicke regime in deriving the spin-spin interaction Hamiltonian, so ground state 
cooling is not necessary [61]. 

We can also obtain effective spin-spin interactions in equation 2 for the case when 
rimO,j <C 6m does not hold, but in this scenario, the phonons in each mode m only 
return to their initial state at particular times l/6m- A favorable situation occurs when 
rjmO,! 6m for m ^ n, that is for all phonon modes except one. In this case, the spin- 
spin interaction takes the form of equation 2 at times tgate = k/6n with k an integer. 
This situation is usually employed when considering quantum gate operations [61]. Note 
that any remaining spin-motion entanglement causes decoherence in the spin state and 
appears as an error in the quantum simulator. 

Using additional sets (two) of Raman helds we can also generate spin-spin couplings 
Tl of the form Hxx + Hyv ~ ® + Sy ® Sy [16, 62]. For that we assume uj is 

tuned close to ± cOm, a situation in which we have spin flip as well as the creation 
of phonons. The expression for is still given by equation 2, but with a re-de£ned 
detuning 6m = uji — (tU'i-i icum)- Note that the different spin-spin couplings, Ej and 
can be independently controlled by the laser parameters of the corresponding Raman 
helds. The corrections to the total Hamiltonian, 

i<j i<j 

arising from the non-commutativity of the different coupling Hamiltonians are oscillatory 
and can be neglected [25, 63, 64]. As experimental complexity increases when considering 
more laser beams, in Appendix D we also discuss a reduced Hamiltonian of the form 
Hxz = Hxx + Hzz, in which the couplings Hyy are not included. Notice that, 
within the perturbation theory on the gauge-invariant manifold, the strong coupling 
Hamiltonians of the two different cases are very similar at lowest orders. We found 
that, in the Hxz case, the underlying physics is almost unaffected when compared to 
the full setup. Alternatively, the simulation may be implemented by a stroboscopic 
sequence of ’’Trotter steps” in which each term in the Hamiltonian is switched on for 
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a short time At consequetively [18]. In this concatenated form, errors dne to non- 
commntativity between the terms in the Hamiltonian scale with the discretized stepsize 
sqnared: 0{AC). This approach also has the advantage that less independent laser 
beams are necessary. 

2.2. Single hexagonal plaguette in an N = 7 ion crystal 

From eqnation (2) it is clear that phonon-mode spectrnm is of central importance to the 
spatial form of the spin-spin interactions, as it determine the Lamb-Dicke parameters 
as well as the mode detnning 6m- In single species ion crystals, in which all ions 
experience the same trapping freqnency, the phonon modes are highly collective, which 
results in long range spin-spin interactions. Here, we discuss how to generate anti¬ 
ferromagnetic spin-spin interactions displaying a hexagonal-plaquette pattern in a 2D 
ion crystal. By plaquette pattern we imply that each of the six spins occupying the 
corners of a hexagon interact with every other spin in the same hexagon with the same 
strength irrespective of their inter-spin separation. First, we discuss how to generate 
it in a minimal setup - a crystal consisting of seven ions - and then extend to a larger 
crystal made of 19 ions. 

The equilibrium conhguration of an ion-crystal with 7 ions in a radially symmetric 
trap is shown in the inset of hgure 2a. Although they form a triangular lattice, we 
picture it as a hexagonal structure with an ion in its centre. The central ion is prepared 
in the atomic state |1), (see hgure 1), which is immune to the Raman helds used for 
spin operations, but experiences an additional one-dimensional pinning optical lattice 
applied normal to the plane of the crystal, i.e. along the 2 ;-axis. The rest of the 
ions occupying the spin states form a hexagon shape. The pinning laser modihes the 
transversal trapping frequency of the central ion to compared to the 

other ions, when the maximum of the lattice potential is focused at its equilibrium 
position, where uol is the harmonic-frequency of the local potential at the ion position 
due to the optical lattice. For a red detuned lattice laser and the ion trapped in 
an antinode, the pinning lattice relaxes the conhnement of the central ion along the 
transversal direction and we see below how it affects the phonon spectrum. In the 
following, we restrict ourselves to pinning-lattice parameters such that the equilibrium 
positions of the ions are not altered once it is switched on. 

The normal-mode spectrum consists of 14 xy (in-plane) and 7 z-(transversal) 
phonon modes, see hgures 3a and 3b. Due to the 2D character of the crystal, the 
in-plane (PM) and the transverse modes (TM) are completely de-coupled. The TMs 
account for the oscillation of ions along the tightly conhned direction, and lie on top of 
the energy spectrum, see hgure 3a. The pinning affects mostly the TMs in two ways: 
(i) it shifts the eigen values and (ii) modihes the mode-structure In particular, the 
frequency ui of the lowest TM (see hgure 2a), which accounts mostly for the oscillation 
of the central ion, moves down to lower values as decreases, see hgure 3c. 

The eigen vector of the lowest TM in the presence of pinning held is shown in 
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Figure 2. (a) The inset shows the equilibirum position of 7 ions in a radially 
symmetric trap with trapping frequencies: 0Jx,y = 27r x 1 MHz and ujz = 2Tr x 3 
MHz. Since the central ion experiences an additional pinning lattice, the 
effective trapping frequency of it along the z axis is reduced to = 27r x 2.7 
MHz. The eigenvector {b\) of the lowest transversal mode of the ion crystal 
is shown in figure 2a. It accounts mostly the oscillation of the central ion, 
and the ions in the outer hexagonal ring oscillates with same amplitude 
and are in-phase, (b) The dimensionless spin-spin couplings Jij/ojQ, where 
uiQ = (^hkj/8M) (uJx/Oj)'^, with (5i = 27r X 10 kHz, the detuning from the lowest 
TM. The slight imperfections from the BFG plaquette interactions arise from 
the off-resonant coupling to higher TMs. For sufficiently small value of Oz, the 
lowest TM mode of the 2D crystal become energetically unstable, i.e. = 0 
and leads to a structural phase transition in which the 2D character of the whole 
crystal has been lost [31]. For the set of parameters in figure 2 it happens when 
Uz < UJzl‘i 2 tt X 1.5 MHz, which restricts the frequency due to the optical 
lattice potential to values ujolI^z < -x/S/d. 


figure 2a. Note that it exhibits a hexagonal character in which all the ions in the 
hexagonal ring oscillate with the same amplitude and in phase. When the Raman beat 
frequency is tuned close to this mode, the resulting spin-spin couplings Ej exhibit a 
hexagonal-plaquette pattern as shown in hgure 2b. There are slight imperfections from 
the ideal hexagonal pattern discussed in the BFG model (equation 1) due to the off- 
resonant couplings between the Raman helds and other TMs. These imperfections can 
be reduced by providing a sufficient gap between the hrst TM and the rest. As shown 
in hgure 3c, this can be done by tuning the frequency of the local harmonic potential 
due to the pinning optical lattice at the central ion position. 

If one considers a Paul trap, the anharmonic coupling between the in-plane micro¬ 
motion and the transverse modes has to be taken into account, which, as recently shown 
[65], only introduces an overall renormalization in the ion positions without effecting 
the normal mode structure. While an overall shift in the transversal mode frequencies 
should also to be taken into account, but it has no effect on our underlying scheme 
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[x 2tt MHz] (jJ^[y. 2 t\ MHz] 

Figure 3. (a) The normal mode spectrum of a trapped 2D ion crystal made 
of 7 ions, in a radially symmetric trap of frequencies Ux^y = 27r x 1 MHz and 
= 27r X 3 MHz. The orange-filled bubbles are the in-plane modes (PM) 
and the blue-filled ones are the transversal modes (TM). (b) The same in the 
presence of the pinning optical lattice. The optical lattice reduced the trapping 
frequency of the central ion along the z axis to a frequency Uz = 2 tt x 2.7 
MHz. The red-filled bubble indicates the lowest TM and the frequency of it is 
shifted to a lower value due to the presence of pinning field (shown by a black 
arrow), (c) The lowest transversal mode vi (filled squares) and the energy 
difference between the lowest two transversal modes <5i^ = 1^2 ~ (filled circles) 
as a function of uiz are shown. For Oz < 0Jzl2 the mode ui gets energetically 
unstable and the 2D nature of the crystal is lost. 


for generating spin-spin interactions. Hence, we can safely neglect those effects in our 
calculations. In the spectrum, there exists a zero-energy in-plane mode due to the 
rotational symmetry of the ion crystal, corresponding to the free rotation of ion crystal 
in the xy plane. In real experiments, any weak stray fields will break this symmtery 
and the crystal will be trapped in a stable configuration. 

2.3. Double hexagonal plaquette in an N = 19 ion crystal 

Having established the basic idea of hexagonal plaquette interactions via an ion crystal 
of seven ions, we expand our discussion to two plaquettes in a larger crystal. Here, 
we consider a crystal of 19 ions in a radially symmetric trap. In order to create two 
plaquettes, we require two pinning standing wave lasers as shown in hgure la and focus 
them such that the maxima of lattice potentials he at the equilibrium positions of the 
two ions populated in state |1). These two ions, as in the previous case, experiences a 
shallower trap along the z axis with respect to the other ions. 

The parameters of the pinning lattices are chosen such that the equilibrium 
positions are not altered. The trapping parameters are taken to he Ux = ojy = 2n x 1 
MHz and = 27r x 3.5 MHz and the effective harmonic frequency at the pinned ions 
along the z axis are = 27r x 2.1 MHz and = 27r x 2.45 MHz respectively 
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for ion 1 and ion 2. The reason for the asymmetry in the pinning lattices parameters 
is discussed below. With these parameters, the two lowest TMs have the plaquette 
character and they comprise mostly of the oscillations of the two pinned ions. If we 
switch off either one of the pinning lattices, we retrieve the same physics discussed in 
the case for one plaquette. As an example, we switch off the pinning lattice on ion 2 
and the ion 1 experiences an effective harmonic trapping frequency of = 27r x 2.1 
MHz. The resulting lowest TM is shown in hgure 4a. It accounts mostly the oscillation 
of ion 1 (not shown in the hgure), and is well localized among the six ions surrounding 
it, forming a hexagonal oscillation pattern. 

The resulting spin-spin interactions, when the Raman helds are detuned close to 
this mode, exhibit the hexagonal-plaquette pattern (see hgure 4b) similar to what we 
have discussed in section 2.2 for the case of 7-ion crystal. This point constitutes one of 
the main results of our paper: even in a larger crystal, by engineering phonon modes we 
can create localized spin-spin interactions, which in our case exhibit a hexagonal pattern. 
Compared to the smaller crystal, here the imperfections are slightly higher. This can 
be understood from the structure of the plaquette eigenmode itself: the amplitude of 
oscillation is not the same for all the six ions in the hexagon (see hgure 4a), which is 
contrary to the N = 7 case. This arises due to the asymmetry in the number of nearest- 
neighbor ions for each ion in the hexagonal ring. A second reason is that the number 
of modes increases linearly with the system size, but in our case the coupling to those 
modes are well suppressed by the large detunings 6m from the corresponding mode. 
Note that if we pin ion 2 instead of ion 1, we create a hexagonal plaquette around the 
second ion. 

Two plaquettes case: To create two plaquettes, we switch on both pinning lattices. 
In this case, the two lowest phonon modes are used to engineer the plaquette interactions. 
If we have identical pinning lattices, the two modes become degenerate. Then, dressing 
them with a single pair of Raman helds may generate an arbitrary superposition of 
phonon states: Cibi ± C 2 b 2 , where bi _2 are the eigen vectors of the two lowest TMs 
and |ci^ 2 p provide us the population of the respective phonon modes. The resulting 
spin-spin interactions do not possess the plaquette character. Hence, we introduce an 
asymmetry between the parameters of the pinning lattices such that the degeneracy is 
lifted as well as the modes are well separated in the spectrum. Then, we use two different 
pairs of Raman helds to induce spin-spin couplings, in which one is detuned near the 
mode bi and the second one is near the mode b 2 with the same detuning 6. Locally, we 
can control the interaction strengths in each plaquette by tuning the ehective coupling 
strengths due to two Raman helds, i.e., 12} and 12}; in the particular example shown 
in hgure 5b we choose 12}/f2} = 1.1 to obtain the same interaction strengths in both 
plaquettes. The equilibrium conhguration of the 2D ion crystal and the eigenvectors of 
the modes bi and b 2 are shown in hgure 5a. In hgure 5a, the ions labelled 13 and 14 
are prepared in the atomic state |1) and pinned by the optical lattices. The resulting 
spin-spin couplings Ej shown in hgure 5b are well localized in each hexagonal plaquette 
and crucially inter-plaquette interactions are negligibly small. Ej constitutes the Ising 
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I I 

Figure 4. (a) The eigen vector for the lowest transversal mode in a 19-ion 
crystal in the presence of a pinning lattice on ion 13. The inset shows the 
pinning scheme in which the ion in the centre of the hexagon is populated in the 
state |1), which experiences an additional optical lattice. The mode accounts 
mostly the oscillation of the central ion (not shown in the plot) and the ions in 
the hexagonal ring. Note that ions outside the hexagon hardly oscillate, (b) The 
dimnensionless spin-spin couplings: Ez jujQ, where ojq = [hkj/8M^ 
for (i = 27r X 10 kHz the detuning from the lowest transversal mode. The 
trapping frequencies are = Wy = 27r x 1 MHz and ojz = x 3.5 MHz and 
the pinned ion experiences an effective shallow trap along the z axis with a 
frequency of = 27r x 2.1 MHz. Note that the interaction pattern is well 
localized in the hexagonal ring. 


part of the Hamiltonian in equation 1. 

The hopping parameter J± is generated by additional Raman fields, again mediated 
through transversal modes. The nature of J± is also vital in determining the ground state 
properties as well as the magnetization dynamics of the emulating spin Hamiltonian. 
As short-range hopping is preferable, e.g. of nearest-neighbor type [8], we identify three 
lowest TMs in the spectrum, labelled as m =1, 2 and 3, see hgure 6a. They lie just 
above the plaquette modes and are also well separated from them, guaranteeing that 
the Ising part is not influenced when addressing with the Raman fields. Note that 
these modes are shared among the two corner-sharing triangles at the centre of the 
crystal, and hence the resulting hopping dynamics may introduce very interesting inter- 
plaquette spin dynamics, for example ring exchange among the four ions {1, 5,15,19} as 
shown schematically in figure 6d, results in charge hopping between the two plaquettes, 
(see section 3 for the consequent magnetization dynamics). The fields are detuned by 
= — 27r X 60 kHz from the m = 1 mode. The calculated couplings are shown in 
hgure 6b and the different hopping channels are schematically shown in hgure 6d. This 
completes our modeling of the microscopic Hamiltonian in an ion crystal for emulating 
the BFG model in equation 1 in a small setup consists of two hexagons. The strategy 
is scalable to larger systems owing to its analogue nature. 
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Figure 5. (a) The eigen vectors of the two lowest transversal modes in a 
trapped 19-ion crystal with trapping frequencies: Wx = Wy = 27r x 1 MHz and 
= 27 r X 3.5 MHz, in the presence of pinning lattices on ions 13 and 14. 
The pinning lattices reduce the effective trapping frequencies of ions 13 and 
14, along the z axis to Lb\ = 2.1 MHz and = 2.45 MHz respectively. The 
inset shows the equilibrium configuration of the corresponding ion crystal. The 
two localized modes: h\ and account mostly the oscillation of ion 13 and 14 
respectively (not shown in the plot), (b) The dimensionless spin-spin couplings: 
J'z with ojQ = (^hk]/8M) for Q'j/O] = 1.1 and 6 = 2tt x 20kHz. 

Note that ion 6 is shared by both the plaquettes. 


3. Many-body Physics and Quantum Magnetism in a double-plaquette 
quantum simulator 

In this section we discuss the many-body physics associated with the above mentioned 
spin systems with hexagonal-plaquette interactions. In particular, we focus on the 
ground state properties and magnetization dynamics. First, we consider the plaquette 
interactions alone, and address the role of the imperfections on the interaction 
patterns at the classical level. Then, we show how gauge-invariance, embodied by the 
magnetization conservation on each plaquette, is affected by the introduction of quantum 
fluctuations, and compare the realistic scenario with the ideal BFG model. Finally, 
we show how interesting many-body dynamics can be observed in a minimal system 
involving two plaquettes, within reach of state-of-the-art experiments with trapped ion 
crystals, by exploring the dynamics of a single charge on top of the gauge-invariant 
background. 

3.1. Ground state properties 

Single plaquette case: In the classical limit {J± = 0) of the Hamiltonian H (equation 1) 
in a single plaquette, the ground state manifold is featured by 20-fold degenerate states. 
Each of them obeys the constraint JZieO = 0) i-®-) f^e total magnetization around 
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Figure 6. (a) The eigen-vectors of three modes lying above the plaquette modes, 
which are used for designing the hopping dynamics, in a trapped = 19 ion 
crystal with trapping frequencies: u)x = u)y = 27r x 1 MHz and = 27r x 3.5 
MHz, in the presence of pinning lattices on ions 13 and 14. The pinning lattices 
reduce the effective trapping frequencies of ions 13 and 14, along the z axis 
to ul = 2.1 MHz and = 2.45 MHz respectively, (b) The dimensionless 
(hopping) spin-spin couplings: Ef for = —27r x 60kHz w.r.t. the mode 
m = 1 shown in (a), (c) shows the ion crystal and (d) shows the different 
hopping channels created by the JV in figure 6b. 


the hexagon is zero or, equivalently, three spins are ”up” and the other three spins are 
’’down”. All such conhgurations are shown in hgure 7. The constraint is akin to the well- 
known ice rule ”2-in 2-out” for the orientation of four dipoles arranged at the corners of 
a tetrahedron in a crystal ice. For the spin-version of the ice setup, so called the ’’spin- 
ice”, the ice rule reads as ”2-up 2-down” [3, 43] and can be mapped into a square lattice 
in which around each vertex four spins are arranged according to ice-rule. There are 
six such possible spin-conhgurations. This six-vertex or ice model required that all the 
four spins around the vertex interact with the same strength, similar to the BFG model 
in which all the spins around the hexagon have the same pairwise interaction strength. 
In our ion-setup, there are imperfections from a perfect plaquette pattern as we have 
discussed above, which lift the 20-fold degeneracy partially. However, as imperfections 
are mostly involving spins within or on adjacent plaquettes, the investigation of a two 
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plaquette system is required to access their effects. 



Figure 7. The 20 different degenerate classical ground state configurations 
which obeys 5'^ = 0 in a single hexagon plaquette. 

Double plaquette case: The number of degenerate ground states increases with the 
system size, and the degree of degeneracy is 200 for two plaquettes with the Ising part 
of the BFG model. In the ion setup, as in the previous case, the imperfections lift this 
degeneracy partially, but the energey levels are closely spaced in the energy spectrum 
as is evident from figure 8a. In our particular example for the ion setup, the classical 
ground states are featured by four degenerate sates and even a very low temperature can 
lead to a classical order-by-disorder phenomenon in which the system collapse into one 
of the ground state configuration. To show the contrast, we compare the above results 
with a nearest-neighbour Ising model in two-plaquettes, with no inter-plaquette spin- 
spin interactions, the ground states are doubly degenerate (antiferromagnetic states) and 
any excited state, resulting from a single spin flip, leads to an energy cost of Jz, see figure 
8a. This implies that our ion setup for the implementation of plaquette interactions is 
indeed a promising candidate for emulating frustrated quantum magnetism and en route 
to quantum spin liquid once scaled up to larger systems. 

In the presence of quantum dynamics the global magnetization, M = SI) is 
preserved, but the magnetization per palquette is not. We define the following operator 
to quantify the magnetization per plaquette. 



where Nq is the number of hexagonal plaquettes (= 2 in our case). The magnitude of 
{G ), the expectation value taken in the ground state of the spin-Hamiltonian, gives us the 
measure of the admixture from states outside the classical ground state manifold obeying 
the plaquette constraint. In practice, this determines the quality of gauge-invariance as 
a function of quantum fluctuations: indeed, G can be viewed as an effective generator of 
a local (gauge) symmetry, which commutes with the effective Hamiltonian in the strong 
coupling Jz ^ J± limit [10]. 

With no hopping (= 0), G = 0, and it increases as increases as shown in figure 
8b. The small values of (G) even for sufficiently large values of J± shows the stability 



{Ei - Eq)/Jz 
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of the low energy manifold with a plaquette constraint against qnantnm fluctnations. 
In both cases, as expected from pertnrbation theory argnments, {G) ~ {J±/Jz)'^ close 
to the classical limit. We note that the only effect of the imperfections is quantitative: 
some gauge variant states not satisfying the plaquette constraint pay a smaller energy 
penalty compare to the ideal BFG case. However, their effect is relatively small even 
for intermediate coupling strengths. 



Figure 8. (a) The lowest 200 classical energy states in a two-plaquette system 
for different spin Hamiltonians with anti-ferromagnetic (AF) interactions. Ei is 
the energy of the state and Eq is the ground state energy. The BFG model 
possesses 200 degenerate states (brown line) obeying the plaquette constraint 
of 3-in 3-out and in the ion setup this degeneracy is lifted by the imperfections, 
but the states are still very closely spaced in the spectrum. For the nearest- 
neighbour (NN) Ising model, there are doubly degenerate AF ground states and 
the next excited states are separated by an energy of J^. (b) The expectation 
value of G in the ground state as a function of the scaled hopping parameter 
J±/ Jz- In the trapped-ions quantum simulator, J±_ is varied by means of Rabi 
frequency (see text). Due to the imperfections in the interactions for the 
trapped ion case, the ratio J±/Jz is defined between their maximum values. 
The small values of (G) even for a sufficiently large values of J_l guarantees us 
the the stability of gauge-invariant low energy manifold. The upper axis for the 
Rabi frequency (corresponding to the J± in the lower axis) is estimated for the 
particular case when Oj = 2tt x 500 kHz, <5^ = 27r x 20 kHz and (5^ = — 27r x 60 
kHz as discussed in section 4. In calculating (G) for the BFG model, we restrict 
the hopping in the five spins located at the intersection of the two plaquettes, 
which makes sense when considering the system in the thermodynamic limit 
and also for the hopping pattern shown in figure 6d for the real situation. 


3.2. Plaguette-magnetization dynamics: deconfined charge excitation 

In this section we examine and compare the dynamics of magnetization in each 
plaquette for an ion setup with plaquette interactions with the BFG model, for different 













CONTENTS 


16 


initially prepared spin configurations. The local plaquette magnetizations are defined 
as Ml = SI and M 2 = respectively for plaquettes 1 and 2. To address 

the stability of gauge constraint against spin hopping, we study the time propagation 
starting from one of the classical degenerate ground state, which obeys 3-in 3-out rule. 
As expected, the hopping dynamics hardly affected the local plaquette magnetizations 
even for large hopping parameter J^. A particular example with J±/Jz = 0.2 is shown 
in figures 9a and 9b for the ion setup and the BFG model, respectively, and they are 
found to exhibit the same behaviour. The dynamics become more interesting when 
we introduce a charge in one of the plaquttes by flipping a single spin in the respective 
plaquette. As an example we consider an initial state in which Mi = 0 and M 2 = — 1 and 
the inter-plaquette spin hopping couples this state to Mi = — 1 and M 2 = 0 and hence 
the state oscillates between them. This has a clear interpretation in gauge theoretical 
language: the ” doped” plaquette represents a charge moving on top of a gauge-invariant 
background. This scenario can be compared to a case where a pair of quasi-particle and 
quasi-hole (fractionalized charges or excitations) is created by a single spin flip in a 
square lattice, which supports an exact incompressible quantum liquid. The particle- 
hole pair doesn’t form a bound pair but is in a deconfined phase [66]. The same is 
also predicted for a classical Ising model in pyrochlore lattice in which the macroscopic 
ground state degeneracy leads to the deconfinement of monopole excitations [67]. This 
is one of the intriguing point of this paper from a frustrated magnetism picture, arising 
with a question on the stability of this deconhned phase, interestingly, which could be 
qualitatively addressed with a small size ion-crystal quantum simulator as we show here. 

The many-body physics associated with the reduced Hamiltonian Hxz mentioned 
in section 2.1 is discussed in Appendix D and found to possesses the same properties as 
that of the full Hamiltonian. 

4. Implementation with a 2D crystal of "‘^Ca^ 

As an example, we consider a crystal of ions. Since this ion does not posses 

a nuclear spin, we employ the Zeeman and metastable electronic states, rather than 
hyperfine ground states. We encode the spins in electronic states of the ion, | ),) = 
= 1/2) and | t) = 1-05/25'"^^ = 3/2) [68]. Inititialization can be performed by 
optical pumping followed by addressed coherent laser pulses that prepare each ion in a 
designated spin state. In this way the ions that will be pinned can also be prepared in 
the state |1) = |S'i/ 2 ,mj = —1/2). The pinning can be done using e.g. a = 532 nm 
retro-reflected laser with a power of ~ 50 mW, focussed to a waist of 2 /im. For an 
estimate, in which we are neglecting the polarizabilities of all transitions except the 
dominant S' —?• P transition around 397 nm, this results in an anti-trapping potential 
of trap-frequencies (0.08,0.08,1.3) 27r MHz in a node of the standing wave. With a Paul 
trap of transverse frequency a;_L = 27r 3 MHz, this results in a local reduction of the trap- 
frequency of ~ 10 % to 27r 2.7 MHz as in the example in section 2.2. The trap-frequency 
in the other two directions is hardly affected (it is reduced by ~ 27r 3 kHz). 
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Figure 9. The time propagation of plaquette magnetizations Mi and M 2 for (a) 
the BFG model and (b) the ion setup with J±/Jz = 0.2 (for the ion setup the 
ratio is taken with the maximum values of and J_l ). The solid lines are for the 
case in which the initial state is prepared in a state in which the Mi = M 2 = 0 
and the magnetization is hardly affected by the hopping dynamics as expected. 
The lines with bubbles are the case for which the initial state is such that one 
of the plaquette carries a charge (in this particular example a negative charge 
in the second plaquette, i.e., M 2 = —1) and the other plaquette has zero charge 
(or zero magnetization, Mi = 0). In time the charge oscillates between the two 
plaquettes indicating a de-confined phase. 


Errors in the realization of the BFG model arise from spontaneous scattering from 
the pinning laser, wich is reduced because of the large detuning and because the ion 
is placed in a node of a wave. We expect a scattering rate of < 1 s“^ as an upper 
bound. The ions are separated by about 5 pm, which is much larger than the waist of 
the pinning laser, such that the ions in the plaquette are not affected by it. 

From an experimental point of view, it is convenient to perform a change of basis to 
exchange az and ax, as spin-spin interactions in the latter require less laser power [69], 
whereas we are interested in the regime Jz > J±- Spin-spin interaction terms involving 
ay can be ommitted in a hrst implementation and this results in non-resonant ® S~^ 
and S~ ^S~ terms that are suppressed by the Gauss law, shown in Appendix D. The Jz 
plaquette spin-spin interaction can now be induced by a Mplmer-Sprensen scheme [61], 
operating on the optical qubit at 729 nm [68, 70]. In particular, a bichromatic light 
held with components around ± Um ± Sm induces Sx ® Sx interactions [69] . Rabi 
frequencies of flj = 27i 500 kHz and a detuning 6m = 2,71 20 kHz result in (plaquette) 
spin-spin interactions of Jz ~ 27r 18.5 kHz. The pinned ion would have a transient 
maximal amplitude of oscillation of around 20 nm -C Xpin, such that the harmonic 
approximation is justihed for the pinning laser. The hdelity of the gate for az is 
calculated semi-classically and its lower bound is about 99.7%, limited by remaining 
spin-motion entanglement in far detuned phonon modes. For details regarding the 
hdelity calculation see Appendix G. To realise J±/Jz ~ 0.2, we require a Rabi frequency 
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for the hopping spin-spin interactions of about Oj ~ 27r 20 kHz with 5^ = —27r 60 kHz 
as discussed above. 

For ions with a nuclear spin, such as ®Be^or the three required long-lived 

states can be encoded in the hyperfine ground states and the spin-spin interactions can 
be induced by Raman lasers, which can reach comparable coupling strengths as for the 
electronic state encoding descibed above [71] . For Yb"*" strong optical pinning has been 
shown recently [55]. 

The readout of the quantum simulation run is performed by detecting state- 
dependent laser induced resonance fluorescence, where the state | J,) will appear bright 
while sites with ions in | '[') remain dark as the state does not scatter photons. Detection 
fidelities exceed 99% in ion trap experiments routinely. Note that the detection can 
easily resolve individual ions, such that the outcome can be immaged on a ccd chip in 
one picture, where populations, and correlations may be readily revealed. This allows 
not only for obtaining the global magnetization but also to determine correlations of 
spins such that the gauge-invariant dynamics can be observed in the plaquettes. 

5. Conclusion and Outlook 

In conclusion, we have shown that using phonon mode shaping it is possible to simulate 
exotic spin-spin interactions in an ion crystal which have far reaching consequences in 
the context of strongly correlated systems, in particular frustrated magnetic systems 
in which low energy manifold is described by an emergent dynamical gauge field. The 
mode-shaping is accomplished by (anti-) pinning standing wave light fields, which modify 
locally the ion trap and affect its motional dynamics. Alternatively, we can selectively 
excite ions to a Rydberg state [57, 58] in which those ions experience a different 
trapping potential compared to the ground state ions due to state-dependent atomic 
polarizabilities, see Appendix E. Additional microwave fields could be employed for 
fine-tuning [59, 60] the trapping field of the ions occupying the Rydberg state. This 
leads to a design of the tranversal mode structure in the planar crystal, similar to that 
with optical pinning forces. 

From a quantum magnetism point of view, we have discussed the ground state and 
magnetization dynamics of the corresponding spin Hamiltonian of the proposed trapped 
ion implementation and the results have been compared with that of the idealized BFG 
model, within a small setup of one and two plaquettes. The results show the excellent 
agreement between the physics in an ion based quantum simulator and that of the 
BFG model. The deconhned dynamics of a single charge excitation in a two plaquette 
system, arising from the ring exchange, is akin to the deconhned monopole excitations 
in a 3D pyrochlore spin-ice, arising from the large macroscopic ground state degeneracy. 
In general, the current studies would open up a completely new aspect of studying 
spin-Hamiltonians using ultra cold ion crystals, in particular once scaled up, will be the 
prime quantum simulator for emulating topological quantum spin liquids or resonating- 
valence-bond states. 
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Appendix A. Atom-light interactions and cr^ G) <Jz interactions 


In this section we discuss briefly the interaction of Raman helds with the spin states 
{| '!'))! t)} cind the resulting spin-spin interactions between the ions. The Hamiltonian 
for N harmonically trapped two level ions interacting with a laser held of frequency cui 
is given by 

H = Ho + Hi, (A.l) 


where 


Ha = 


hw 


N 


m=3N 


It 






(A.2) 

2=1 m=l 

where is the Pauli spin-1/2 matrix and the operator am {^In) annihilates (creates) a 
phonon in the mode with an eigen value Um (see Appendix B). In the interaction 
picture, the atom-light interaction Hamiltonian [11, 72, 73] reads as. 






(A.3) 


where Q] = —dS^{2 G M is the effective coupling strength for ion with d, the dipole 
moment for the transition and is the electric held strength at the ion position, 
dm = dmoKn^i ^^re the Lamb-Dicke parameters with = ■\/h/2Mz/m and = {&m} 
is the phonon eigen vector, and 0/’s are additional phases in each ion. The operator k*, 
a general 2x2 matrix, is determined by the polarization of the electric held and the 
atomic states. Note that the summation of m in the exponential term in equation A.3 
is only over the N transversal 2 ;-modes and we replaced the position operators of the 
ions in terms of normal mode operators while writing the Hamiltonian Hj [11]. 

In the weak coupling limit: Qj and when uj (an implementation in 

this limit is shown in hgure lb) and in the Lamb-Dicke regime: {{dm + dm)^) ■C 1 

where we keep only the terms upto hrst order in the Lamb-Dicke parameters r]l^, then 
taking the rotating wave approximation (RWA) in which the terms with are 

neglected, we hnally get the atom-light interaction Hamiltonian as 

N 

Ht = E + hx., (A.4) 

i m=l 

where 5m = — i^m- From the form of equation A.4 it is clear that there is no spin 

flip involved which is justihed by the limit uji -C wi-i-. In general, equation A.4 describes 
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quantum gates based on an effective laser field being interacting with many ions. The 
electric field from the laser beams gives rise to a Stark shift for each internal states, 
and a corresponding electric-dipole force on each ion. These state dependent forces are 
typically used in trapped-ion quantum computation [61, 74]. The time evolution of the 
system is then obtained using the Magnus expansion [75] for the time evolution operator 


U{t, 0) = exp 





1 



Jo 


(A.5) 


where H{t) is the time-dependent Hamiltonian. In our case, with sufficiently large 
detuning from the motional sideband [16] and assuming the same phase for all ions, the 
long-term time evolution is dominated by a term linear in time, t and is of the form 

| 15 ]with 


N 


4' = E 




m=l 


(A.6) 


This realizes an Ising spin-Hamiltonian of the form = Ylij ^ with spin-spin 
couplings Ej between and ions. Note that the spin-1/2 operators (So) defined in 
equation 1 are related to the Pauli spin-1/2 matrices (cTq,) through Sa = 


Appendix B. Ion crystal: Equilibrium positions and vibrational spectrum 


In a 3D setup, the potential experienced by the ions is 


i=l 


STren r 


(B.l) 

r^'ji &r\ 1— r’-1 

where |rj — r^l = ^{xi — Xj)^ + {Vi — DjY + {zi — ZjY- The first term is the harmonic 
confinement due to external fields and the second term accounts for the Coulomb 
interaction between the positively charged ions. We choose the length scale as = 
(7^e^)/(47reoMa;^) and in the dimensionless form the potential reads as 


= E 4 h? + + w t 


j-r. (B.2) 

r,'— rJ 

1=1 * ij=i^i^j 

where = Uz/oOx, \ = ojy/ux ■ In the following we absorb in the position co¬ 
ordinates such that Xi, Ui and Zi are dimensionless. The equilibrium positions are 
calculated by solving the equations 


dV _dV _dV 
dxi dyt dzi 

for each ions. This leads to solving the coupled algebraic equations of the form 

N 


\2 ^ I ~ _ n 






(B.3) 


(B.4) 
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where a = {x,y,z} and-Ro(hj) = {xi — XjY + iyi — yjY + ^Zi — ZjY. Once the equilibrium 
positions are obtained, we can calculate the eigen vectors and eigen values of 
the phonon modes by the exact diagonalization of the Hessian matrix constructed out 
of the following second order derivatives: 


for m = n 


ioT m ^ n 


d^V 


N r 


dcX'indcXm 

d^V 

0(Xndf3rn 

d^V 




3 ( 0 ;^ - araf _ 1 


N 




3 ( 0 :* Oim){/3i Pm) 




Ro{i, mp'R 

3(cin CUm) 


dcx^dcXm 

d^V 
d(Xnd Pm 

with a,P E {x,y,z} and p. 


+ 


Ro{n,mYR RQ{n,m)^R 

3(Q:n ^m){,Pn Pm 

R^in^mpO 


(B.5) 


(B.6) 


(B.7) 

(B.8) 


Appendix C. Semi-classical estimation of population in the motional states 
and gate fidelity 


The phonon spectrum, 


Hm — hVm ^ ^ ^\rflm 


(C.l) 


is described as a collection of harmonic oscillators with frequencies Vm- As described in 
the text, the quantum emulation of the BFG Hamiltonian relies on the phonon-mode 
shaping by optically anti-pinning few ions in the crystal. Our focus will be on low lying 
transverse modes, which gives us a hexagonal plaquette pattern when coupled dressed 
by Raman laser helds. Any discrepancy from the exact BFG model in the emulation 
results from the population of other modes which is not of our prime interest. We semi- 
classically investigate the population of these higher modes, and provide an estimation 
for the fidelity of the corresponding gate operations. The problem is equivalent to an 
undamped, driven harmonic oscillator in which the driving force is a sinusoidal one with 
frequency cu/, that of the laser held. The governing equation is 


(Rz 2 


hQiko 
M 


cos{ujt), 


where ko = rjm/qmo with 
z'{0) = 0, we get the solution 

hfljko 


h 


2Ml'rr 


(G.2) 

For the initial conditions ;2(0) = 0 and 


z{t) = 


M {Rm - ^ 1 ) 


COS OJit — cos Umt] . 


(G.3) 
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Figure Cl. The percentage of fidelity loss as a function of the displacement 
zlqmo- To keep the fidelity loss with in 1%, the displacement zjqmo required 
to be roughly <0.25. 


The classical, (maximum) amplitude of oscillatiou, A ~ J/iOjfco related 

to the quautum mechauical spectrum of Harmouical oscillator by equatiug the classical 
aud quautum euergies, aud we get 




QmO J 


(C.4) 


after we left out the zero poiut oscillatiou euergy. This relatiou would provide us a 
rough estimatiou for the uumber of phouous iu the system for a particular mode m. lu 
additiou, the overlap betweeu the iuitial phouou state aud the maximally displaced state 
would give us the hdelity loss arisiug from this particular mode. We require this loss 
would be well below 1%. Takiug the Gaussiau grouud state solutiou of the harmouic 
oscillator, we get the overlap as a fuuctiou of the displacemeut zq as 


Om = exp(-2;2/8g^o). 


(C.5) 


The hdelity loss Cm = 1 — Om due to a giveu mode m, as a fuuctiou of the amplitude 
of oscillatiou is showu iu Fig. Cl, aud oue cau see that the larger the displacemeuts, 
the lesser the hdelity is. Iu order to keep the hdelity loss with iu 1%, the displacemeut 
z/qmo required to be <0.30. The uet hdelity loss is theu obtaiued by the summiug over 
that of all modes. 


Appendix D. Gauss’s law and plaquette magnetization dynamics for the 
reduced Hamiltonian, Hxz 

In this section we consider a reduced Hamiltonian of the form 

= + (D.l) 

i<j i<j 

Experimentally, it is advantageous to have the reduced Hamiltonian over the full one 
since less number of independent lasers are involved. The coupling parameters Ej and 
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JlNz 


Figure Dl. (a) Gauss’s law for the reduced Hamiltonian Hxz compared to 
the full Hamiltonian H in an ion-setup for the same parameters as in figure 
8. (b)The time propagation of plaquette magnetizations Mi and M 2 for the 
reduced Hamiltonian with Jz = 0.4 (the ratio is taken with the maximum 
values of Jz and J_l). The initial state is prepared such that Mi = 0 and 
M 2 = -1. 


El are still the same which we discussed in the main part of the text. Contrary to the 
full Hamiltonian, here it involves not only flip-flop terms but also flip-flip and flop-flop, 
as seen from the expansion, SxSx = {S+S- + S-S+ + S-S- + S'+S'+)/4. The results 
are shown in hgure Dl and the resulting Gauss’s law is compared with that of the full 
Hamiltonian H. The hopping of charge excitation is also demonstrated, revealing the 
high flexibility of trapped-ion quantum simulators. 

Appendix E. Rydberg excitations in an ion crystal: atomic polarizibility 
and modified trap frequency 

Instead of using the pinning lattices on ions in state |1) (see main text), we can 
as well excite or weakly dress those ions to a Rydberg state. Due to the state 
dependent atomic polarizabilities, the Rydberg excited ions experience different trapping 
frequencies compared to those occupying the low-lying atomic states. This modihes the 
phonon spectrum of the ion crystal, identical to the situation in which the pinning 
lattices are present. Recently, it has been shown [59, 60] that ions excited to the nP- 
Rydberg state experience an additional radial potential, Vaijc) ^ —e^a^VnpE, where 
the polarisability VnP ~ —0.25 x rC atomic units, with n, the principal quantum 
number of the Rydberg state. Effectively, an ion in this Rydberg state experiences a 
tighter conhnement compared to an ion occupying a low-lying state, but with additional 
microwave helds, by coupling to a nearby Rydberg state, say n'S, we can freely tune 
the trapping frequency of the Rydberg excited ion. It can be increased, decreased or 
even made equal to that of ions occupying the qubit states | t) and | J,). The dressed 
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states |±) = N±{C±\nP) + Iri'S')) polarizabilities read as 

V± = Nl{CiVnp + Vn's) (E.l) 

with Vn's > 0, the parameter C± depends on the microwave parameters such as 
detunings and Rabi frequencies and N± are the normalization constants. More details 
on the microwave-Rydberg approach can be found in [60]. V± are easily controlled via 
microwave held paremeters, and hence the local trapping of Rydberg ions. The only 
limit in this approach is set by the life time of Rydberg states for gate operations and 
it can be augmented by going either to a higher Rydberg level or by weakly dressing to 
the Rydberg state. 
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